Experimental implementation of heat-bath algorithmic cooling using solid-state 

nuclear magnètic resonance 



Jonathan BaughQ Osama Moussa, Colm A. Ryan, Ashwin NayakQ and Raymond LafLammc[| 
Institute for Quantum Computing, University of Waterloo, Waterioo, Ontarïo N2L 3G1 

(Dated: Nov. 24, 2005) 

The counter-intuitive properties of quantum mechanics have the potential to revolutionize in- 
formation processing by enabling efficient algorithms with no known classical counterparts Q, ■ 
Harnessing this power requires developing a set of building blocks £g|, one of which is a method 
to initialize the set of quantum bits (qubits) to a known state. Additionally, fresh ancillary qubits 
must be available during the course of computation to achieve fault tolerance QIEISQ· ^ n an y 
physical system used to implement quantum computation, one must therefore be able to selectively 
and dynamically remove entropy from the part of the system that is to be mapped to qubits. One 
such method is an " open-system" cooling protocol in which a subset of qubits can be brought 
into contact with an external large heat-capacity system. Theoretical efforts 0, 0, lïo| have led to 
an implementation-independent cooling procedure, namely heat-bath algorithmic cooling (HBAC). 
These efforts have culminated with the proposal of an optimal algorithm, the partner-pairing algo- 
rithm (PPA), which was used to compute the physical limits of HBAC [ïï|. We report here the first 
experimental realization of multi-step cooling of a quantum system via HBAC. The experiment was 
carried out using nuclear magnètic resonance (NMR) of a solid-state ensemble three-qubit system. 
It demonstrates the repeated repolarization of a particular qubit to an effective spin-bath tempera- 
ture and alternating logical operations within the three-qubit subspace to ultimately cool a second 
qubit below this temperature. Demonstration of the control necessary for these operations is an 
important milestone in the control of solid-state NMR qubits and toward fault-tolerant quantum 
computing. 



I. INTRODUCTION 

NMR-based ensemble quantum information processing 
(QIP) devices have provided excellent testbeds for con- 
trolling non-trivial numbers of qubits 0, 0, 0, 0] . A 
solid-state NMR QIP architecture builds on this success 
by incorporating the esscntial features of the liquid-state 
devices while offering the potential to reach unit polar- 

ization and thus control more qubits mm. In this ar- 
chitecture, the abundant nuclear spins with polarization 
P form a large heat-capacity spin-bath that can be either 
coupled to, or decoupled from, a dilute, embedded ensem- 
ble of spin-labelled isotopomers that comprise the qubit 
register. Bulk spin-cooling procedures such as dynamic 
nuclear polarization are well known and capable of rcach- 
ing polarizations near unity 0, 0] . This architecture is 
one realization within a large class of possible solid-state 
QIP systems in which coherently controlled qubits can 
be brought into contact with an external system that be- 
haves as a heat bath. The principies and methods applicd 
in solid-state NMR QIP will therefore apply to many 
other systems. An additional motivation is development 
of control techniques that future quantum devices will 
utilize. For this experiment, we develop a novel technique 
to implement the controlled qubit-bath interaction, and 
also report the first application of strongly-modulating 
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pulses 18| to solid-state NMR for high-fidelity, coherent 
qubit control. 



II. THREE-QUBIT MALONIC ACID SYSTEM 

The three-qubit quantum information processor used 
here is formed by the three spin-1/2 13 C nuclei of iso- 
topically labelled malonic acid molècules, occupying a 
dilute fraction of lattice sites in an othcrwise unlabeled 
single-crystal of malonic acid (unlabeled, with the ex- 
ception of naturally occurring 13 C isotopes at the rate 
of 1.1%). The concentration of labelled molècules was 
3.2%. Malonic acid also contains abundant spin-1/2 1 H 
nuclei, which comprise the heat-bath. Figure^shows the 
1 H-decoupled, 13 C-NMR spectrum for the crystal (and 
crystal orientation) used in this work. The spectrum 
shows the NMR absorption peaks of both the qubit spins 
(quartets) and natural abundance 13 C spins (singlets), 
the latter being inconsequential for QIP purposes. The 
table in Fig.^lists the parameters of the ensemble qubit 
Hamiltonian obtained from fitting the spectrum, and also 
includcs couplings involving the methylene protons cal- 
culated for this crystal orientation from the known crys- 
tal structure [19|. Experiments were performed at room 
temperature at a static magnètic field strength of 7.1 T, 
where the thermal 1 H polarization is Pji — 2.4 x 10 -5 . 
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over 13 C, 1 H nuclei, respectively, and cr| is the /3-axis 



FIG. 1: Characteristics of the dilute 3- 13 C malonic acid spin 
system. (below) 1 H-decoupled, 13 C spectrum near the [010] 
orientation with respect to the static magnètic field. The blue- 
dashed line is the experimental NMR absorption spectrum, 
and the solid-red line is a fit. Multiplet assignments are indi- 
cated by the labels Ci, C2 and C m . The central peaks in each 
multiplet correspond to natural abundance 13 C in the sam- 
ple, which are inconseqüent ial for QIP purposes. The peak 
height difïerences in the 3- 13 C molecule peaks indicate the 
strong coupling regime, i.e. the 13 C- 13 C intramolecular dipo- 
lar couplings are significant compared to the relative chemical 
shifts. (above) Table showing the 13 C rotating-frame Hamil- 
tonian parameters (chemical shifts along diagonal; dipolar 
coupling strengths off-diagonal; all vàlues in kHz) obtained 
from the spectral fit. It also includes calculated dipolar cou- 
plings invol ving the methylene protons based on the atòmic 
coordinates pj|| and the crystal orientation obtained from the 
spectral fit. 



III. REFRESH OPERATION 

In this orientation, the methylene carbon C m has a 
dipolar coupling of 19 kHz to H m \ of the methylene 1 H 
pair, whereas no other 13 C- 1 H dipolar coupling in the 
system is larger than 2 kHz. Therefore, a spin-exchange 
Hamiltonian of the form 



E 

j£C,k£H 



(1) 



that couples the two nuclear species will generate dynam- 
ics dominated by the large C m -H m i coupling at short 
times (the Djk are 13 C- 1 H dipolar couplings, indices run 



Pauli operator for spin a). 
coupling Hamiltonian, 
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we applied a multiple-pulse 'time-suspension' sequcncc 
[20| synchronously to both 13 C and X H spins to create 
the effective spin-exchange Hamiltonian (in the toggling 
frame), to lowest order in the Magnus expansion of the 
average Hamiltonian [2Ï|. Application of the sequence 
for the C m -H ml exchange period r = 4xl g kHz ~ 40/is 
results in an approximate swap gate (state exchange) 
between the C m and H m \ spins. With an initial bulk 
l· l·l polarization Pjj, this procedure yields a selective 
dynamic transfer of polarization P' = tjPh to C m , where 
< i) < 1 and ideally 1 77 1 = 1 . We define the effective 
spin-bath temperature to be that which corresponds to 
the experimentally obtained P' under this procedure, 
and refer to this transfer as a refresh operation. We 
obtained P' ~ 0.83Pf/ experimentally, and found that 
repeated refresh operations showed no loss in efRciency 
given at least a 6 ms delay for 1 H- 1 H equilibration. 
However, we observed a decay of Pr as a function of the 
number of repetitions, due to accumulated control errors, 
which lcad to an identical loss in the refresh polarization. 



IV. IMPLEMENTATION OF 
PARTNER-PAIRING ALGORITHM 

The experiment consists of the first six operations of 
the partner-pairing algorithm (PPA) on three qubits: 
three refresh operations, and three permutation gates 
that operate on the qubit register. This is described 
in the quantum circuit diagram of Fig. [3 During the 
register operations, the l· l·l polarization is first rotated 
into the transverse plane, and then 'spin-locked' by a 
strong, phase-matched RF field that both preserves the 
bulk l· l·l polarization and decouples the 1 H- 13 C dipolar 
interactions. Since 1 H- 1 H dipolar interactions are 
merely scaled by a factor —1/2 under spin-locking, H m i 
is allowed to equilibrate with the bulk l· l·l nuclei via spin 
diffusion. This occurs on a timescale longer than the 
transverse dephasing time (T2(H m ) ~ 100/is), but much 
shorter than the spin-lattice relaxation time (Tf 1 ~ 50s) 
of H m \. Hence, H m i plays the role of the fast-relaxing 
qubit described in the protocol of Schulman et al. [ïï| . 
The first two register operations are swap gates: the 
third is a three-bit compression (3BC) gate 0, 0, E3 
that boosts the polarization of the first qubit, Ci, at 
the expense of the polarizations of the other two qubits. 
Ideally, the protocol builds a uniform polarization on all 
three qubits corresponding to the bath polarization (first 
five steps), then selectively transfers as much entropy as 
possible from the first qubit to the other two (last step). 
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FIG. 2: Schematic quantum circuit diagram of the imple- 
mented protocol. Time fiows from left to right. The three-bit 
compression (3BC) gate is shown here decomposed as control- 
not gates and a control-control-not (Toffoli) gate. The gate 
sequence corresponds to the first six steps of the partner- 
pairing algorithm |TÏ|| on three qubits. The input state is a 
collective polarization of the bulk l· l·l spins. The refresh oper- 
ation is approximately 40/ís in duration, whereas the register 
operations are between 0.7 and 1.3 ms in duration. Thermal 
contact takes place during 1 H spin-locking pulses that begin 
just prior to the register operations, and extend an additional 
12 ms after each operation. H m i can be thought of as an 
additional 'special purpose' qubit in this experiment; despite 
non-selective 1 H control (due to bulk hydrogenation), the re- 
fresh and thermal contact operations could be performed us- 
ing collective 1 H control. Thus, H m i serves as a fast-relaxing 
'qubit' and the bulk 1 H-bath as a large heat-capacity thermal 
bath. 



The last step (3BC) leads to a polarization boost by a 
factor of 3/2 on the first qubit. Subsequently, the heated 
qubits can be re-cooled to the spin-bath temperature, 
and the compression step repeated, iteratively, until 
the asymptotic value of the first-bit polarization is 
reached. This limiting polarization depends only on the 
number of qubits and the bath polarization and is 
ideally P(Ci) = 2P' for three qubits (for n qubits it is 
2«-2p/ m regime P' « 2~ n , and 1.0 in the regime 
P' » 2~ n (refs. [Ulli). The first six steps carried 
out here should yield a polarization of 1.5P' on Ci, 
assuming ideal operations. 

The control operations performed herein are quan- 
tum control operations: state-independent unitary 
rotations in the Hilbert space. However, it should be 
noted that the HBAC gates are all permutations that 
map computational basis states to other computational 
basis states. Therefore, gate fidelities were measured 
with respect to correlation with these known states, 



rather than the manifold of genèric quantum states. 
We took advantage of this property to further opti- 
mize the control parameters of the 13 C gates (register 
operations) for the state-specific transformations of 
the protocol. These operations were carried out using 
numerically optimized control sequences referred to 
as strongly-modulating pulses fisj ]. Such pulses drive 
the system strongly at all times, such that the average 
RF amplitude is comparable to, or greater than, the 
magnitude of the intcrnal Hamiltonian. This allows 
inhomogeneities in the ensemble qubit Hamiltonian to 
be efficiently refocused, so that ensemble coherence is 
better maintained throughout the gate operations. 

In this set of experiments, the 13 C qubit spins are 
initialized to infinite temperature (a preceding broad- 
band 13 C ir/2 excitation pulse is followed by a dephasing 
period in which 1 H dipolar fields effectively dephase the 
13 C polarization). Following the fifth step, polarizations 
(in units of P') of 0.88, 0.83 and 0.76 (±0.03) are built 
up on Ci, Ci and C m , respectively. The final 3BC 
operation yields P(d)/P' = 1.22 ± 0.03, a boost of 48% 
compared to the average polarization (0.82) following 
step five. Despite control imperfections that effectively 
heat the qubits at each step, we are able to cool the 
Ci qubit ensemble well below the effective 1 H spin-bath 
temperature. 



V. RESULTS 

The results are summarized in Fig. |3| in (a) are 
shown the spectral intensities corresponding to 13 C spin 
polarizations following each of the six steps, and in (b) 
the integrated intensities are graphed in comparison 
with the ideal vàlues. We note that the overall fidelity 
of the experiment, F = 1.22/1.50 = 0.81, implies an 
error per step of 3.7%. This error rate is only about a 
factor of two larger than the average error per two-qubit 
gate obtained in a benchmark liquid-state NMR QIP 
experiment [ï^. Furthermore, the state-correlation 
fidelity of the 3BC gate over the polarizations on all 
three qubits is 0.96 ± 0.03. From Fig. EO», it can be 
seen that the fidelity of the refresh operation drops off 
roughly quadratically in the number of steps; this is 
consistent with the loss of bulk 1 H polarization due to 
pulse imperfections both in the multiplc-pulse refresh 
operations and in the spin-locking sequence. Since the 
broadband pulses have been optimized for flip-angle in 
these sequences, we suspect that the remaining errors 
are mainly due to switching transients that oceur in the 
tuned RF circuitry of the NMR probehead, and to a 
lesser extent off-resonance and finite pulse-width effects 
that modify the average Hamiltonian pçj. Similar 
effects lead to imperfect fidelity of the 13 C control. 
With suitable improvements to the resonant circuit 
response and by incorporating numerical optimization 
of the multiple-pulse refresh operations, we expect 
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that several iterations of the protocol could be carried 
out and that the limiting polarization of 2P' could be 
approached in this system. The same methodologies 
should also be applicable in larger qubit systems with 
similar architecture. For a 6-qubit system using the 
PPA, a bath polarization P > 0.2 would be sufhcient, 
in principle, to reach a pure state on one qubit |22| . 
Such bulk nuclear polarizations are well within reach 
via well-known dynamic nuclear polarization techniques 
[Ï7|: for example, unpaired electron spins at defects 
(g = 2) in a field of 3.4 T and at temperature 4.2 K are 
polarized to 0.5. 

This work demonstrates that solid-state NMR QIP 
devices could be used to implcment active error cor- 
rection. Given a bath polarization near unity, the 
refresh operation implemented here would constitute the 
dynamic resetting of a chosen qubit. This would allow 
a new NMR-based testbed for the ideas of quantum 
error correction and for controlled open-system quantum 
dynamics in the regime of high state purity and up to 
~ 20 qubits. 



VI. METHODS 

NMR experiments were carried out at room temper- 
ature on a Bruker Avance solid-state spectrometer op- 
crating at a field of 7.1 T, and home-built dual channcl 
RF probe. The sample coil had an inner diameter of 
3 mm, and the employed tt/2 broadband pulse lengths 
were 1.25/us and 0.75/is for 13 C and 1 H, respectively. The 
sample was a 4.5 x 1 x 1 mm 3 single crystal of malonic 
acid grown from aqueous solution with a 3.2% molecu- 
lar fraction of 3— 13 C labelled molècules. Spectra were 
obtained by signal averaging for 80 scans. The proton 
spin-lattice relaxation time was Tf 1 = 50s, so the delay 
between scans was set to 62~i = 300s. Design of the 
strongly-modulating 13 C pulses followed very closely the 
methodology described in ref. , and penalty functions 
were adjusted to favor average RF amplitudes compara- 
ble to or greater than the magnitude of the 13 C rotating- 
frame Hamiltonian. These pulses were optimized and 
simulated over a 5-point distribution of RF amplitudc 
corresponding to the measured distribution over the spin 
ensemble (cr = 6.2% in RF amplitude). The 'time- 
suspension' sequence applied synchronously to 13 C and 
1 H was a 12- pulse subsequence of the Cory 48-pulse se- 
quence [2ÏJ] ■ The delays between pulses were adjusted so 
that the total length of the sequence was 40/is. 1 H spin- 
locking/decoupling was carried out at RF amplitude of 



250 kHz. The spectra in Fig.[3]were obtained by applying 
a 7r/2 broadband pulse to read out the spin polarizations. 
The readout pulses were preceded by a 1.5 ms delay in 
which 1 H decoupling was on but any off-diagonal terms in 
the 13 C spin density matrix would significantly dephase 
(T 2 * ~ 2ms). The absolute value of the refresh polariza- 
tion P' was determined by comparing the initial refresh 
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FIG. 3: Experimental results in terms of C spectra and 
their integrated intensities. (a) Readout spectra obtained fol- 
lowing each of the six steps in the protocol. The integrated 
peak intensities for each multiplet correspond to the ensem- 
ble spin polarizations. The natural abundance C m signal that 
appears at each refresh step (adding to the intensity of the 
central peaks) should be ignored; we are only interested in 
the part of the signal arising from the 3- 13 C qubit molècules, 
which can be seen clearly in the Ci and C2 spectral regions, 
(b) Bars indicate ideal qubit polarizations at each step; exper- 
imental vàlues obtained from integration of the above spectra 
are shown as shaded bands, whose thickness indicates exper- 
imental uncertainty. 



polarization on C rn and the thermal equilibrium 13 C po- 
larization Pq measured in a separate experiment. These 
yield the ratio of P' to Pc, and P' to Pr using the fact 
that P H = 3.98P C - 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation [2] Quant. Inform. Process. 3, 1 (2004), special issue. 

and Quantum Information (Cambridge University Press, [3] D. P. DiVincenzo, Fort. der Phys. 48, 771 (2000), quant- 
Cambridge, UK, 2000). ph/0002077. 



5 



[4] E. Knill, R. Laflamme, and W. H. Zurek, Science 279, 
342 (1998). 

[5] A. Y. Kitaev, in Quantum Cornmunicatïon, Computing 
and Measurement (Plenum, New York, USA, 1997), pp. 
181-188. 

[6] D. Aharonov and M. Ben-Or, in Proc. 29th. Ann. ACM 
Symp. on Theory of Computing (1997), longer version 
quant-ph/9906129 quant-ph/9611025. 

[7] J. Preskill, Proc. R. Soc. Lond. A 454, 385 (1998), quant- 
ph/9705031. 

[8] L. Schulman and U. Vazirani, Proc. of the 31th Annual 
ACM Symposium on Theory of Computing pp. 322-329 
(1999). 

[9] P. O. Boykin, T. Mor, V. Roychowdhury, F. Vatan, and 
R. Vrijen, Proc. Natl. Acad. Sci. USA 99, 3388 (2002), 
quant-ph/0106093. 

[10] J. M. Fernandez, S. Lloyd, T. Mor, and V. Roychowd- 
hury, Int. J. Quant. Inform. 2, 461 (2004). 

[11] L. Schulman, T. Mor, and Y. Weinstein, Phys. Rev. Lett. 
94, 120501 (2005). 

[12] E. Knill, R. Laflamme, R. Martínez, and C.-H. Tseng, 
Nature 404, 368 (2000), quant-ph/9908051, URL 
http : //www. ari v . org/ abs/quant-ph/9908051 

[13] N. Gershenfeld and I. L. Chuang, 275, 350 (1997). 

[14] D. G. Cory, M. D. Price, and T. F. Havel, Physica D 
120, 82 (1998), quant-ph/9709001. 

[15] D. G. Cory, R. Laflamme, E. Knill, L. Viola, T. F. Havel, 
N.Boulant, G. Boutis, E. Fortunato, S. Lloyd, R. Mar- 



tinez, et al., Fort. der Phys. special issue, Experimental 
Proposals for Quantum Computation 48 (2000), quant- 
ph/0004104. 

[16] G. M. Leskowitz, R. A. Olsen, N.Ghaderi, and L. J. 

Mueller, J. Chem. Phys. 119, 1643 (2003). 
[17] A. Abragam and M. Goldman, Nuclear Magnetism: Or- 

der and Disorder (Oxford University Press, Oxford, Eng- 

land, 1982). 

[18] E. M. Fortunato, M. A. Pravia, N. Boulant, G. Tekle- 

mariam, T. F. Havel, and D. G. Cory, J. Chem. Phys. 

116, 7599 (2002). 
[19] N. R. Jagannathan, S. S. Rajan, and E. Subramanian, J. 

Chem. Cryst. 24, 75 (1994). 
[20] D. G. Cory, J. B. Miller, and A. N. Garroway, J. Mag. 

Res. 90, 205 (1990). 
[21] U. Haeberlen, High Resolution NMR in Sòlids: Selective 

Averagmg (Acadèmic Press, New York, USA, 1976). 
[22] O. Moussa, On heat-bath algorithmic cooling and its im- 

plementation in solid-state NMR (2005), MSc. thesis, 

ht t p : / / www .iqc.ca/ omoussa/ work /thesis . 
We gratefully acknowledge: D. G. Cory, T. F. Havel, and C. 
Ramanathan for discussions and use of NMR simulation code; 
W. P. Power, M. Ditty and N. J. Taylor for facility use and 
experimental assistance; ARDA and NSERC for support. O. 
M. acknowledges the Ontario Ministry of Training, Colleges 
and Universities for support. 

Correspondence and requests for material should be sent to 
J. B. (e-mail: baugh@iqc.ca). 



